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Using numerical simulations, this work demonstrates a concept called “back-end 
ignition” for lighting-off and pre-heating a catalytic monolith in a power generation system. 
In this concept, a downstream heat source (e.g. a flame) or resistive heating in the 
downstream portion of the monolith initiates a localized catalytic reaction which 
subsequently propagates upstream and heats the entire monolith. The simulations used a 
transient numerical model of a single catalytic channel which characterizes the behavior of 
the entire monolith. The model treats both the gas and solid phases and includes detailed 
homogeneous and heterogeneous reactions. An important parameter in the model for back- 
end ignition is upstream heat conduction along the solid. The simulations used both dry and 
wet CO chemistry as a model fuel for the proof-of-eoncept calculations; the presence of 
water vapor can trigger homogenous reactions, provided that gas-phase temperatures are 
adequately high and there is sufficient fuel remaining after surface reactions. With 
sufficiently high inlet equivalence ratio, back-end ignition occurs using the thermophysical 
properties of both a ceramic and metal monolith (coated with platinum in both cases), with 
the heat-up times significantly faster for the metal monolith. For lower equivalence ratios, 
back-end ignition occurs without upstream propagation. Once light-off and propagation 
occur, the inlet equivalence ratio could be reduced significantly while still maintaining an 
ignited monolith as demonstrated by calculations using complete monolith heating. 


Nomenclature 

a* = ratio of catalytic to geometric surface area 

= channel open flow area (m 2 ) 

= solid cross-sectional area (m 2 ) 

= C p k Y k ; mean gas specific heat at constant pressure (J / kg /K) 

= species k gas specific heat at constant pressure (J / kg /K) 

= solid specific heat (J / kg /K) 

= species k diffusion coefficient into mixture (m 2 /s) 

= enthalpy of species k (J/kg) 
hok = convection mass transfer coefficient for specie k in bulk mixture (m/s) 

hr = internal convection heat transfer coefficient (W / m 2 K) 
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species index 

K g +K S ; total number of species (gas + surface) 
forward rate constant for the /th reaction (units m, mole, s) 
number of gas-phase species 

reverse rate constant for the /th reaction (units m, mole, s) 

thermal conductivity of solid (J/m/s/K) 

number of surface species 

mass flow rate down the channel (kg/s) 

number of gas-phase reactions 

number of surface reactions 

Upstream face Nusselt number for monolith 

rate of progress of the / th gas reaction (mole / m 3 /s) or 

rate of progress of the / th surface reaction (mole / m 2 /s) 

universal gas-constant (8.314 J / mol / K) 

net production rate of specie k due to all surface reactions (mole / m 2 /s) 

circumferential length of channel cross section = nd (m) or; 
surface area per unit length (m) 
time (s) 

time integration between successive gas-phase integrations (s) 
bulk gas temperature (K) 

temperature of species k entering or leaving control volume 
(if s k > 0 then T[ =T S otherwise T[ =T) 

solid temperature (K) 

average (bulk) velocity of flow (m/s) 

molecular weight of specie k (g/mole) 

-l 

; mean molecular weight of the mixture (g/mole) 

axial coordinate from channel inlet (m) 
molar concentration of specie k (mole / m 3 ) 
mass fraction of specie k in the bulk mixture 
mass fraction of specie k at the wall 
site fraction of surface specie k 




Greek Letters 

F = surface density of site on solid (2.7063 x 10" 5 mol / m 2 for platinum) 

v ki = net stoichiometric coefficient of the k th specie in the /th reaction; 

v’ki = stoichiometric coefficient of the k th reactant specie in the /th reaction 

v” ki = stoichiometric coefficient of the k th product specie in the /th reaction 

p = mass density of gas (kg/m 3 ) 

p s = mass density of solid (kg/m 3 ) 

CT k = site occupancy number - number of sites that species k occupies on surface 

4> = equivalence ratio 

6) k = gas-phase production rate of specie k due to all gas reactions (mole/m 3 /s) 


I. Introduction 

C atalytic combustion has long been recognized for its potential to reduce temperatures of power generation 
devices in both large and small scale applications including microcombustors. Such power generation devices 
often use a channel geometry (or multiple channels such as in a monolith structure) coated with a catalyst through 
which flows a pre -mixed combustible gas mixture. The catalyst must be at an elevated temperature to initiate a 
reaction. One challenge for such devices is the startup, or light-off, from a cold initial condition. We propose a new 
concept for the light-off of a catalytic monolith utilizing a “back-end” ignition scheme. In this concept, only a small 
portion of the downstream catalyst section is heated (via a gaseous flame stabilized downstream of the monolith or 
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Joule heating of the solid) initiating a localized 
catalytic reaction (see Fig. 1). The catalytic 
reaction releases sufficient heat which is 
conducted upstream via the solid causing the 
reaction zone to propagate, ultimately heating 
the entire monolith. Recently, we have 
predicted such a phenomenon using a numerical 
model and verified the predictions with an 
experiment in an isolated non-adiabatic single 
tube made of pure platinum. 1 " 3 * 

In this work, we utilize the same numerical 
model, which is of a transient catalytic channel, 
but extend the results to a monolith 
configuration. The fuel used for the simulation 
is premixed carbon monoxide in air. CO is an 
excellent model gas for which detailed chemical mechanisms (both homogeneous and heterogeneous) are readily 
available in the literature but are not so complex as to make transient computations intractable. The inlet flow 
conditions are consistent with small scale power generation requirements while the inlet equivalence ratio is kept in 
a range to prevent the surface temperature from exceeding material property limits. 

This paper examines the characteristics of both ceramic and metallic monolith substrates which differ 
significantly in thermophysical properties and hence have different operating requirements for startup. In 
commercial monoliths, the substrate material typically is coated with a catalytic washcoat, which depending on the 
amount of catalyst added (i.e. catalytic loading), may enhance the catalytic activity by increasing the available 
surface sites per unit length of the solid. The model accounts for catalytic loading by defining a parameter a* which, 
in effect, increases the geometric surface area per unit length available for catalysis. This parameter, which can 
enhance the amount of catalytic conversion, can account for the variety of catalyst loading conditions encountered in 
commercial monoliths. 


Back-end Ignition 


Steady state 
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Homogeneous 
flames to promote 
catalytic reactions 


Stabilized catalytic 
reactions on surface 



Catalytic reactions propagate upstream 

Figure 1. Back-end ignition concept for monolith reactor 
using downstream homogeneous flame. 


II. Model Description 


A. Overview 

A complete description of the model is presented in Ref. 3 and only a top-level overview is presented here. The 
model is of a single catalytic channel unit (gas plus solid section) and represents an interior channel of a monolith 
that can be approximated as adiabatic. The single channel model is used to characterize the overall performance of 
the monolith. The model considers both the gas and solid phase of the channel. A premixed gas of fuel and air 
enters the channel at a prescribed velocity, temperature, and pressure. The model assumes that the channel is made 
of a solid substrate material (either ceramic or metal) which is coated with a thin catalytic material (washcoat). The 
combustion process considers both detailed gas phase and surface reactions. There is negligible pressure drop along 
the channel (i.e. no momentum equation solved) and the gas obeys the ideal gas law as shown in Eq. (1). 

Equations (l)-(4) describe the gas-phase model which is quasi-steady relative to the transient solid, Eqs. (5)-(7), 
owing to the significantly longer heat-up times of the solid (~sec) compared with the residence times of the gas in 
the channel (~ msec) 3 4 . Equation (2) states that there is no net consumption of mass from the gas phase. In the gas- 
phase equations, heat and mass diffusion in the axial direction are neglected because the Peclet number based on 
typical gas velocities is much greater than unity. The model captures the heat transfer perpendicular to the flow 
using a heat transfer coefficient, h T . In Eq. (3), if the reaction rate of species k due to catalytic reactions, j , is 

greater than 0 (gas-species k desorbing from the surface) then the enthalpy of species k (h k ) is evaluated at the 
surface temperature (T' = T s ). For j < 0 (gas-species k adsorbing to the surface), h k is evaluated at the bulk gas 

temperature (T' = T). 

There are K g gas-phase species and the model uses two values of gas-phase species mass-fraction, a bulk flow 
value (Y k ) and a value adjacent to the catalytic surface but still in the gas phase (Y kw ). Lateral mass transfer 
coefficients come from the heat transfer coefficients using the analogy of heat and mass transfer 5 . The value of Y kw 
comes from Eq. (5) which balances lateral mass transfer and the adsorption or desorption of species on the catalytic 
surface. Transport properties are mass averaged at each axial location at the bulk temperature. 

Equation (6) shows the transient solid-phase energy equation which includes heat transfer to and from the gas 
inside the tube (term 2), heat generation terms due to catalytic reactions (terms 3 and 4), internal heat generation 
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(term 5), and axial heat conduction along the solid (term 6). The enthalpy of absorbing and desorbing gas species 
are evaluated just as described previously for the gas phase. For both gas and surface species, the enthalpy comes 
from data provided with the chemical mechanisms. Because the solid is thermally thin, the surface catalytic 
reactions (those involving surface species only) are modeled as surface heat generation in the solid. A prescribed 
volumetric heat generation (term 5) simulates Joule or external heating of the solid. The solid density, p s , heat 
capacity, C s , and thermal conductivity, k s , are all constant in the model. The model accounts for vaiying surface 
coverage of adsorbed species. There are K s total surface species including vacant surface sites. Equation (7) 
accounts for the time variation of surface species and assumes that the total number of active surface sites, T, 
remains constant. The value for T (=2.7063xl0' 5 mol / m 2 for platinum) comes from data provided with the 
heterogeneous mechanism. 

B. Equations 

1. Gas phase 


Equation of State: 



Overall mass-conservation: m = pnA = constant 


(2) 


Energy conservation: 


r)T Kg Kg 

puACp — + A-£ a> k W k h k + s k W k [h k (T) - h k (T )] + h T S[r -T s ] = 0 

OX k = 1 k = 1 

r \T ifi A .<0 
[T s if s K > 0 


(3) 


Conservation for species k: 


puA d -^ = s k W k S + co k W k A 

OX 


(4) 


2. Surface and Solid-Phase Equations 


Flux-matching at surface: 

ph» k (Y k -Y kw ) = -s k W k 


(5) 

Energy conservation: 

Ps^s A s = h T S(r — T s ) 

ot 

d 2 T 

+ <jge„ A S +kA S _ 0 
OX 

-S%s k W t h k (T') + S%s t W t h k (T s ) 

k = 1 k= 1 

(6) 

Rate of change of surface 
site fractions: 

dZ. 1 . ,, , v , 

— - = -s k (k= l...K s ) 

dt r 


(7) 

3. Chemistry’ Expressions 

n„ 

1=1 


(8) 
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The formulae for evaluating the gas-phase, co , and surface, t , reactions are shown in Eqs. (8)-(10) where q; is 

the rate of the i th reaction and v is the stoichiometric coefficient. Formulas for the forward and reverse reaction rate 
constants may be found in the literature (for example, see Ref. 10). To account for an area enhancement due to the 
use of a thin catalytic washcoat, the geometric surface area per unit length, S, is multiplied by a surface area 
adjustment factor, a*, for terms dealing with surface reactions but not terms dealing with gas-to-surface transport. 
The parameter a* can also be thought of as the ratio of effective catalytic surface area to the geometric area. While 
a* is thought of here as a surface area enhancement, it is mathematically equivalent to modifying each kinetic pre- 
exponential term by the factor a*. In this study, a* is constant and set to a value of 20 which was previously shown 
to accelerate the kinetics sufficiently so that the reactions, once ignited, were mass-transfer limited. 3 

C. Solution Procedure 

In general, the dependent variables (T. T s , Y k , Y kw , and Z k ) as well as all the gas property values, transport 
coefficients, and reaction rate terms are functions of axial position and time. Since the gas phase is quasi-steady, it 
responds instantly to changes on the solid surface. Thus, the gas phase must be solved simultaneously with the 
changing solid phase and surface variables. This fact is especially important for cases which are mass-diffusion 
limited. The model separately integrates the gas-phase equations, Eqs. (3)-(4) which are ODEs, along the length of 
the reactor at specific time intervals. For this integration in x, the spatial distributions of the solid and surface 
parameters (T s , Z k , Y kw , and S k ) are fixed at the particular instant in time. The integration takes the necessary 

spatial step to handle stiff regions (i.e. ignition) interpolating the surface parameters, as necessary. 

Equations (5)-(7) represent the surface / solid-phase equations and are directly coupled to the gas phase via the 
variables Y k and T. The method of lines transforms Eq. (6) from a PDE into a system of ODEs. The solid is 
discretized into many control volumes; Equations (5)-(7) are written separately for each solid control volume thus 
forming a large system of ordinary differential-algebraic equations (DAE) for the solid. 

The inclusion of detailed gas-phase and surface chemistry introduces numerical difficulty via mathematical 
stiffness. Stiffness in DAE systems essentially means that the time (or spatial) step required to solve the equations is 
much smaller than that required to obtain an accurate solution. Thus, stiff DAE systems require special numerical 
routines to solve them efficiently. The model uses the publicly available code DASPK 5 6 , which is specifically 
designed for large systems of stiff equations, to perform the integrations. 

The basic solution algorithm is to integrate the solid / surface equations, Eqs. (5)-(7), a finite amount forward in 
time, At, assuming that the axially varying gas-phase values are constant parameters during the interval At. The 
code, upon reaching At, integrates the gas-phase equations, Eqs. (3)-(4), along the length of the channel keeping the 
solid / surface values (T s , Z k , Y kw , and s k ) fixed. This method of solution essentially lags the gas phase by the 

amount At during the time integration. For the calculations presented in this paper, At was set to approximately 0.03 
sec. (or about l/100 th of a characteristic solid time scale based on interior gas-solid heat transfer time scale of ~3 sec. 
for the conditions in this paper) and provided time-accurate results. More details on the integration process are 
presented elsewhere. 3 The series of spatial and temporal integrations continued to steady state which was based on 
observing that the temperature profiles were no longer changing to within the numerical accuracy requested of the 
solution. 

D. Chemistry 

The fuel considered is carbon monoxide (CO) and includes both dry and wet versions - the wet version includes 
hydrogen chemistry which is used here to allow for the presence of water vapor in the feed stream. In the dry 
mechanism, the model tracks the following species: O, CO, 0 2 , C0 2 , N 2 , Pt(s), CO(s), C0 2 (s), C(s), and O(s). The 
“(s)” denotes species adsorbed to the surface while the symbol Pt(s) denotes vacant sites on the platinum surface. 
With wet CO chemistry, the model also includes the following species: Fl 2 , H 2 0, FI 2 0 2 , FI, OH, HCO, H0 2 , H(s), 
H 2 0(s), and OH(s). 

The gas-phase CO mechanism is from Davis et al. 7 For dry CO chemistry, the gas phase utilizes only steps 8, 
21, and 23 of the mechanism - the remaining reactions deal with hydrogen chemistry and are neglected for dry CO 
computations. For calculations using wet CO chemistry, the full mechanism is utilized. 

The heterogeneous reactions, which account for species adsorption, surface reactions, and desorption, come from 
a subset of a CH 4 /0 2 on platinum mechanism proposed by Deutschmann et al. (the website version is used) 8 
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Calculations for dry CO use steps 4, 5, 6, 7, 15, 16, 17, 18, 23, and 24 of the mechanism. For wet CO, all of the 
reaction steps are used except for steps 19-22 which involve the adsorption and subsequent decomposition of CFL|. 

E. Thermophysical Properties and Geometric Parameters 

For the gas phase, the thermodynamic properties (C P and h K ) come from data provided with the mechanism. 
Transport properties come from the Gordon-McBride database or kinetic theory for pure species while the mixture 
uses appropriate mixing rules (see Ref. 3 for more details). For the solid, this paper presents calculation for 2 
representative monolith reactors: one made of cordierite and the other from a metal alloy. Table 1 shows the solid- 
phase property variables used for the computations as well as the geometric parameters of the monolith (which are 
representative for the given material). The property values come from Appendix B of Ref. 9. While there is some 
variation in thermophysical properties of the solid with temperature (especially Cs and k s ), the computations used 
the constant values shown in Table 1 for these proof-of-concept calculations. 


Table 1. Thermophysieal and geometric parameters 


Material 

Ceramic 

Cordierite 

Metal 
70-80% Fe 
10-20% Cr 
5-10% A1 

p s , Density (kg/m 3 ) 

2510 

7400 

C s , Fleat Capacity (J/kg/K) 

1110 

460 

k s . Thermal Conductivity (W/m/K) 

2 

19.4 

d, Inner Diameter (mm) 

1 

1 

A/A t , Open Area Fraction per Monolith Cell 

67% 

88% 

L, Monolith Channel Length (mm) 

60 

60 


F. Initial and Boundary Conditions 

The gas-phase equations, which are first-order, require only a single upstream boundary condition. The inlet 
conditions are representative of a flow downstream of a compressor in a small scale gas-power turbine. Specifically, 
the gas-phase boundary conditions upstream of the monolith are a temperature of 300K (in actual applications the 
temperature may be somewhat higher downstream of the compressor), a pressure of 3 atm., and a fuel equivalence 
ratio ranging from <J)=0. 1 to 0.4 using mixtures of CO and air. The channel flow velocities were either 5, 10, or 20 
m/s and, because of the small channel diameter, the flow remained laminar. For the cases using wet CO chemistry, 
water vapor (1% by volume) is added to the inlet feed. The solid energy equation, Eq. (6), requires 2 boundary 
conditions due to the axial conduction term. There is heat loss to the inlet face of the monolith to the cold inlet gas. 
An estimate of the heat loss 3 comes from stagnation point theory which yields a Nu w = 33, 47, and 66 for an 
upstream velocity of 5, 10, and 20 m/s, respectively, for the ceramic monolith geometry. For the metal monolith 
geometry, the upstream Nuu is 62, 88, and 125, respectively. The outlet of the solid channel is adiabatic. During the 
back-end ignition phase, an externally applied heat source is imposed in the last 10% of the channel length for a 

specified time duration. This is implemented through the internal energy generation term, q , in Eq. (6). 

The transient equations, Eqs. (6)-(7), require initial conditions for surface temperature and site fraction 
distribution along the length of the channel. The computations all began with a cold initial condition. The initial 
surface temperature is 300K while all surface sites are occupied initially by CO (i.e. surface site fraction Z C o(si % ) 
owing to the higher sticking probability of CO relative to 0 2 , the 2 major reacting species in the flow. The initial 

condition for the remaining unknown, Y kw , comes from an algebraic solution of Eq. (5). 

In the numerical simulations, the internal heat-generation term, q , provides the energy to raise the monolith 

temperature to initiate catalytic reactions. The q n is either applied uniformly along the entire monolith or only 

across a portion of the monolith nearest to the outlet (termed back-end ignition). In both cases, heating by q can 

simulate resistive heating in the solid. For the case of back-end ignition, q can simulate heating by an external 

homogeneous flame that was established previously as shown in Fig. 1. In the simulations, q is turned off after a 

given interval of time. In a real monolith, however, the homogeneous flame can continue to burn downstream of the 
monolith (in the after-bed space) provided that there is sufficient reactants (and high-enough gas temperature) even 
after catalytic reactions are established in the monolith." Table 2 shows the parameters used in both the complete 
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reactor heating and the back-end ignition scheme. For back-end ignition, the ceramic required more energy 
deposited to reach an ignition temperature because of the ceramic’s higher heat capacity. 


Table 2. Ignition Parameters for Numerical Simulation 


Ignition Strategy 


Duration (s) 

Heating Entire Monolith: ceramic & metal 

30 

1 

Back-end Ignition (last 10% of length): Ceramic 

3 

1 

Back-end Ignition (last 10% of length): Metal 

3 

0.4 


ITT. Results 

The results include both steady-state and transient data. The steady-state results include the fuel conversion as 
well as temperatures (both solid and gas) at the exit of the channel. For select cases, the species and temperature 
profiles along the channel are presented. The transient data presented is the propagation speed of the catalytic 
reactions front which can be used to gauge the preheat time of the channel. 

A. Steady-state data 

Figure 2 shows the fuel conversion (defined as the mass of gas-phase fuel reacted to the mass of gas-phase fuel 
fed into the reactor) at the outlet of the channel as a function of inlet equivalence ratio for the ceramic monolith. 
The calculations using the metal monolith produced very similar results in steady state and are not presented. In the 
figure, three separate computations are shown for each velocity: (1) complete monolith heating (labeled “All”) with 
dry CO chemistry, (2) back-end monolith heating (labeled “Back”) with dry CO chemistry, and (3) complete 
monolith heating using wet CO chemistry. In all cases, the internal heat generation was turned off after 1 second. 
As seen in Fig. 2, the model did not light-off (i.e. no steady conversion after removal of the heat source) for either 
heating profile at a <(> < 0. 1. At cf) = 0.11, heating the entire monolith allowed light-off to occur and to be sustained 
resulting in fuel conversion at steady state. With heating the entire monolith, an increasing 4> produced an increasing 
conversion with the largest increase (15% to 25%) occurring from 4> = 0. 1 1 to 0.15. Between <|> = 0. 15 and 0.4 there 
was only a modest 10-15% increase in total conversion. 

For back-end ignition (dotted lines in Fig. 2), a higher (|) of 0.15 was required to see sustainable conversion in the 
channel. From <j) = 0.15 to 0.23, the light-off position occurred at the channel exit and remained near the exit 
yielding only a small total conversion. For (|) > 0.23, the light-off position initially occurred at the outlet but 
subsequently propagated toward the inlet. In steady state, the resulting channel conversion for the cases which 
propagated (i.e. (|) > 0.23) matched that seen with complete channel heating. 

In Fig. 2, the exit conversion decreased with increasing velocity. For instance, the conversion at (|) = 0.4 was 
98%, 81%, and 51% for u = 5, 10, and 20 m/s. As the channel velocity increases from 5 to 20 m/s, the gas-phase 
residence time decreases from 12 ms to 4 ms, respectively. For a given channel diameter, the mass-transport time 
from the gas to the surface is more or less fixed (~2 ms for the present channel, see Ref. 3 and 4 for more 
discussion). Since the reactions are essentially mass-transfer limited in these cases, the faster velocity cases do not 
have enough residence time in the channel to completely react. 

Figure 3 shows the solid-phase temperature (corresponding to the same calculations in Fig. 2) at the outlet of the 
monolith. At (j)=0.1 and below, there is no light-off and the steady-state temperature is the ambient value of 300K. 
Above <|> =0.1, the temperature at the end of the monolith increases with increasing <|>. For complete monolith 
heating, the monolith’s exit temperature ranges from -650K at <|> =0.11 to - 1 63 OK at <|> =0.4. Using back-end 
ignition, the channel exit temperature is -100K lower compared to temperature achieved with complete monolith 
heating between <|> =0.15 and 0.23. For (j) > 0.23, the catalytic reaction front propagated upstream and the steady- 
state solutions are the same irrespective of the heating profile. The temperature limits shown in Fig. 3 are taken 
from Ref. 9 and are exceeded for continuous operation with <|> > 0.35 for this particular monolith configuration. 

Figure 4 shows the gas-phase temperature at the outlet of the monolith. The trends in this figure are the same as 
those seen in Fig. 3 except that the peak temperatures are somewhat less than the solid-phase since the gas is being 
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Figure 2. Steady-state fuel conversion through the ceramic monolith as a function of inlet equivalence 
ratio for a CO / air inlet feed. For the wet CO cases, water-vapor (1% by volume) is added to the inlet 
feed. Two different types of heating profiles are used in the simulation: complete monolith heating 
(“All”) and back-end heating of the monolith (“Back”). 




Figure 4. Steady-state gas-phase temperature at the exit of the ceramic monolith. 

8 


American Institute of Aeronautics and Astronautics 


Complete Monolith Heating Back-End Ignition 



Figure 5. Steady-state temperature and species profiles along the length of the ceramic monolith channel 
after complete and back-end monolith heating. The inlet conditions are (|)=0.2 at 10 m/s. 

heated by the solid. For ([> > 0.35 and u=10 m/s, the computations using wet CO show a higher exit temperature (as 
well as a higher conversion as seen back in Fig. 2) compared with calculations using dry CO. In the wet CO case, 
gas-phase reactions have been initiated due to the presence of water vapor which provides hydroxyl (OH) radicals at 
the elevated temperatures. This particular condition is somewhat unique among the computations in that there is 
enough fuel remaining when gas temperatures become sufficiently high to initiate gas-phase reactions. For the 
u=20m/s, the residence time is too short to heat the gas-phase sufficiently for homogeneous reactions to occur. For 
the u=5m/s cases, the gas-phase temperatures near the outlet are sufficiently high at (|) ss 0.27 and greater. In these 
cases, however, catalytic reactions have consumed most of the fuel, leaving little fuel for homogeneous reactions. 

Figure 5 shows the steady-state temperature (top most graphs) and species profiles (middle and bottom graphs) 
along the length of the ceramic monolith. The inlet flow conditions are 4> = 0.2 at 10 m/s. At these conditions, two 
different steady states occurred depending on the monolith heating profile. For complete monolith heating (left 
column), significant catalytic reactions begin at the inlet as is evident by the high solid temperature and low CO 
mass-fraction (at the wall) exactly at the inlet. The solid-phase inlet temperature profile shows that there is upstream 
heat-loss given the non-zero temperature gradient at the inlet. The surface reactions are mass-transfer limited as is 
suggested by the large (almost 3 orders of magnitude) difference between the gas-phase bulk and wall mass- 
fractions for CO (reactions can be considered mass-transfer limited if there is a large difference between the bulk 
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and wall values 9 ). Regarding the surface species, O(s) 
is the major surface species in steady state with 
complete monolith heating. 

For the case of back-end ignition (right column of 
Fig. 5), the catalytic reaction front anchors itself by the 
channel outlet near where the ignition energy was 
applied (see region labeled Ignition Zone in Fig. 5). In 
this case (4> = 0.2 at 10 m/s), the catalytic reaction did 
not propagate upstream. The solid-phase temperature 
near the outlet was over 900K but the gas was only 
slightly heated at the outlet (~ 350K) due to the short 
contact time with the hot surface. The wall-mass 
fractions of CO suggested that once catalytic reactions 
initiated, they became mass-transfer limited across only 
a short distance. For this case, the major surface 
species along the channel was CO(s) until the light-off 
point (x » 5.5 cm). After light-off, the major surface 
species was O(s). These data suggest that up until 
light-off, the reaction is limited by the lack of O(s) on 
the surface. After light-off, the reactions appear to be 
limited by the low concentration of CO(s) on the 
surface which is likely a consequence of the slow 
transport rate (relative to the potential reactions rate) of 
gaseous CO to the surface. 

For the (j)=0.2 case shown in Fig. 5, the gas-phase 
temperature rises from 300K to 900K across the length 
of the channel was insufficient to trigger gas-phase reactions. Previously, Figs. 2 and 4 showed that for inlet 
conditions of (|)=0.4 and u=10 m/s, there was a higher gas-phase exit temperature and greater CO conversion when 
comparing the dry and wet CO. For the wet CO cases, water vapor in the amount of 1% (by volume) was added to 
the inlet gas. Figure 6 shows the steady-state mass-fractions of the gas-phase species for the <j)=0.4 and u=10 m/s 
calculations. The species common to both dry and wet CO calculations are shown in color. The species present 
only in the wet CO calculations are shown using hatched lines. For the dry CO calculations, only the major stable 
species are present in large quantities and there is very little O which indicates that no significant gas-phase 
reactions are occurring. The formation of CO 2 is due almost exclusively to surface reactions. For the wet CO 
calculations, the data in Fig. 6 show the presence of significant gas-phase radicals (especially O and OH) caused by 
the addition of water vapor in the inlet feed. The presence of water vapor (which introduces the OH radical) lowers 
the temperature required to initiate gas-phase reactions for CO / 0 2 systems, as is well established in the combustion 
literature (for example, see Ref. 10). Comparing the bulk CO mass-fraction profiles (blue lines), it is clear from Fig. 
6 that gas-phase reactions are causing further reduction of bulk CO (i.e. more conversion) in the wet CO calculations 
compared to the dry CO. 

B. Transient Results 

To this point, the results have shown that different steady states can occur depending on initial heating profile. 
With regard to back-end ignition, the cases of practical interest are those in which the catalytic reaction propagates 
upstream and heats the entire monolith. For these cases, a parameter of interest to the designer of catalytic monolith 
reactors is the heat-up time of the solid. Figure 7 shows the transient evolution of the solid-temperature (as well as 
the final steady state of both the solid and gas-phase temperature) for both a ceramic (left graph) and metal (right 
graph) monolith. The data in these figures correspond to inlet conditions of 4>=0.4 and u=10 m/s and show both dry 
and wet CO computations using the back-end ignition scheme. For the ceramic monolith, the ignition term (q ) 

was on for 1 second. For the metal monolith, the igniter duration was 0.4 seconds and was reduced to avoid 
unnecessarily high solid temperatures during the ignition process. Figure 7 shows that, just after removal of the heat 
generation term, the solid temperature was -720K and 850K for the ceramic and metal monolith, respectively. It 
takes the reaction front about 78 seconds to propagate about midway across the channel for the case of a ceramic 
monolith and requires about 171 seconds to achieve steady-states conditions. For the metal monolith, it takes the 



Distance down the channel (cm) 

Figure 6. Steady-state mass-fractions comparing 
the dry and wet CO calculations at inlet conditions 
of (|i=0.4 and u=10 m/s. The colored lines represent 
species in both the dry and wet CO computations. 
The hatched lines are onlv nart of wet CO. 
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temperature front only 15 seconds to propagate half-way across the channel and only about 34 seconds to achieve 
steady state. 



Distance down the channel (cm) 

Figure 7. Transient evolution of the solid-phase temperature along the channel for both a ceramic (left 
graph) and metallic substrate (right graph). The computations are for an inlet condition of <j>=0.4 and 
u=10 m/s and include both dry and wet CO cases. 

The propagation velocity of the temperature front was nearly constant along the channel especially away from 
the ends of the channel. The temperature front was defined as the first location (starting from x=0) where the solid 
was 650K or greater. The propagation velocity, V P , comes from a linear curve fit of the position versus time across 
the region from 25% to 75% of the propagation distance. The propagation distance was defined as the region 
between the location of the temperature front when the heat-generation was turned off (about x=5.5 cm in Fig. 7) 
and x=0. 



Equivalence Ratio 

Figure 8. Propagation velocity of the temperature front as a function of inlet equivalence ratio and velocity 
during back-end ignition. 


Figure 8 shows the propagation velocity of the solid-phase temperature front versus equivalence ratio for both 
the ceramic and metal monoliths. The propagation speed for the metal monolith is significantly higher than for the 
ceramic. For the ceramic monolith at a 10 m/s inlet velocity, the propagation speed of the reaction front ranged 
from V P = 0.07 mm/s at (j)=0.25 to Vp= 0.35 mm/s at (j)=0.4 (with only slight variations with flow velocity). For the 
metal monolith at 10 m/s inlet velocity, the propagation speeds increased to V P = 0.54 mm/s at <(> = 0.25 and Vp = 
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1.83 mm/s at (J> = 0.4 (with slightly greater variations with flow velocity as shown in the figure). The propagation 
speed was roughly 7 times faster in the metal compared with the ceramic for the <|> = 0.25 cases. This ratio is close 
to the ratio of metal to ceramic thermal diffusivities (~8). For an equivalence ratio of 4> = 0.4, the propagation speed 
was about 5.25 times faster relative to the ceramic. The propagation speed ratios varied only slightly with inlet 
velocity. Interestingly, the equivalence ratio for which propagation occurred is not that different between the metal 
and ceramic monolith. The temperature front propagated at approximately 4>=0.24 for the ceramic and only slightly 
leaner, (j)=0.225, for the metal. 


IV. Discussion 

The data from the numerical simulations suggest that back-end ignition of a catalytic monolith for power 
generation devices is a feasible concept. With only a few watts of power (per channel) delivered to the back end of 
the monolith (say from a flame or resistive heating) to initiate localized catalytic reactions, the chemical energy 
released during the subsequent upstream propagation of the catalytic reactions can provide the remaining energy to 
heat the entire monolith. It is important for the reaction front to propagate upstream so that the solid temperature is 
hot along the majority of the channel. This allows the gas sufficient time to heat as it passes through the monolith 
which is the objective of the monolith (i.e. combustor) for power generation devices. 

There is a reasonable range of equivalence ratios, about (j)=0.25 to 0.35, for which back-end ignition and 
propagation occurred for the monolith configuration tested in the simulations. For this reaction, the lower 4> 
corresponds to the minimum equivalence ratio for which propagation occurred while the highest 4> corresponds to 
the material temperature limits. Since the data in this paper is specific to the monolith parameters used in the 
simulations, further work which can generalize when back-end ignition occurs is recommended. There is likely 
some combination of heat generation (due to catalytic reactions), heat-transfer to the flowing gas, and upstream heat 
conduction along the solid which allows the reaction front to propagate. Such a correlation would be valuable to the 
designer of a catalytic monolith reactor using the back-end ignition concept. 

For a practical combustor using back-end ignition, the operating sequence can be more elaborate than tested in 
the simulations. During the ignition sequence, for instance, the inlet flow can be initially throttled back (e.g. the 5 
m/s data) as well as operated under richer conditions. This may be important for establishing a downstream flame 
(such as shown schematically in Fig. 1) to avoid ignitability problems and blow-off. Once the catalytic reactions are 
initiated, the mixture can be made leaner to avoid exceeding any material property limits. In this situation, the 
downstream flame would eventually extinguish due to both the leaner inlet conditions as well as due to reactant 
consumption by the catalytic reactions. Once ignited and preheated, the operating range in terms of inlet equivalence 
ratio for the monolith reactor increases (e.g. down to almost (j)=0.1 as demonstrated by the fully heated monolith 
cases). 

Under certain conditions ((j)=0.4 and 10 m/s as an example for the monolith configuration in the present work), 
gas-phase reactions can be initiated within the monolith. This can be an advantage since no further catalyst is 
required downstream of the onset of homogeneous reactions. The reaction can go to completion in the gas phase 11 
reducing the length of the monolith which is advantageous from a cost and perhaps a design standpoint. 
Furthermore, limiting the length of the catalytic monolith prevents further heating of the solid if homogeneous 
reactions occur near the outlet. As a possible disadvantage, there is a risk of initiating homogeneous reaction too 
early within the monolith (if inlet conditions are mismanaged) potentially causing excessive temperatures 
downstream. Also, in many applications, the main emphasis of using a catalytic monolith is to avoid the high- 
temperatures of gas-phase reactions for both pollution control and material property limits in the power generation 
device. Fortunately, the data presented in this paper suggests that the onset of homogenous reactions occurs in only 
limited conditions and can be controlled with proper management of inlet conditions and monolith dimensions. We 
note that catalytically induced gas-phase reactions downstream of the monolith have been previously demonstrated 
experimentally. 11 


V. Conclusion 

Using numerical simulations, this work demonstrates a concept called back-end ignition for lighting-off and 
heating a catalytic monolith in a power generation system. This concept uses a heat source (e.g. a flame) or resistive 
heating in the downstream portion of the monolith initiating a localized catalytic reaction which subsequently 
propagates upstream, heating the entire monolith. The simulations used a transient numerical model of a single 
catalytic channel. This channel represents an interior channel of the monolith which is approximated as adiabatic 
and is used to characterize the behavior of the entire monolith. The model treats both the gas and solid phases and 
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includes detailed homogeneous and heterogeneous reactions. The simulations used chemical mechanisms of both 
dry and wet CO and air mixtures as a model fuel for the proof-of-concept calculations. The solution procedure takes 
advantage of the slow solid time constants relative to the gas phase allowing the use of a quasi-steady approximation 
for the latter. An important parameter in the model for back-end ignition is upstream heat-conduction along the 
solid. 

Back-end ignition occurred using both a ceramic and metal monolith, with the heat-up times significantly faster 
for the metal monolith. In some cases, back-end ignition occurred without upstream propagation. For an inlet 
velocity of 5 m/s at 3 atm., back-end ignition occurred at 4> > 0.15 but propagation of the reaction front required a 
richer inlet equivalence ratio (<|> > 0.24). For both the higher inlet velocities tests (10 and 20 m/s), slightly richer 
conditions (4> > 0.25) were required for upstream catalytic propagation. The propagation speed of the reaction front 
along the metal monolith ranged from about 0.5 to 1.8 mm/s (with some variation with inlet velocity) for (|)=0.25 and 
0.4, respectively. For the ceramic monolith, the propagation velocities were significantly slower ranging from -0.05 
to 0.36 mm/s for the same equivalence ratio range. Once light-off and propagation occurred, the inlet equivalence 
ratio could be reduced significantly (e.g. (|)«0.11) and still maintain an ignited monolith as demonstrated by 
calculations using complete monolith heating. 

Under certain conditions (cj)=0.4 and 10 m/s as an example for the monolith configuration in the present work), 
gas-phase reactions were predicted to occur at the back end of the monolith. In these computations, significantly 
more CO conversion occurred with the presence of water vapor (1% by volume) in the inlet feed. The water vapor 
introduced the OF! radical which subsequently lowered the temperature required to initiate gas-phase reactions with 
CO and air. Besides having OF! radicals and high enough gas-phase temperatures, sufficient fuel must remain in the 
gas phase to initiate gas-phase reactions. At long residence (slower inlet velocities), catalytic reactions alone can 
consume most of the fuel in the reactor. 
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